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Confinement of spin-orbit induced Dirac states in quantum point contacts 
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The quantum transmission problem for a particle moving in a quantum point contact in the 
presence of a Rashba spin-orbit interaction and applied magnetic field Is solved semiclasslcally. A 
strong Rashba interaction and parallel magnetic field form emergent Dirac states at the center of the 
constriction, leading to the appearance of resonances which carry spin current and become bound at 
high magnetic fields. These states can be controlled in situ by modulation of external electric and 
magnetic fields, and can be used to turn the channel into a spin pump which operates at zero bias. 
It is shown that this effect is currently experimentally accessible in p-type quantum point contacts. 
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One dimensional (ID) spin-orbit coupled systems have 
recently attracted significant interest in the context of 
quantum information and spintronics, playing a key role 
in the search for emergent Majorana fermionsir— and 
the generation of spin polarized current^^— . Interest in 
these systems has sparked several theoretical studies of 
their modified conductance propertiesiSr— . In this work 
I present a semiclassical solution to the scattering prob¬ 
lem for a quantum point contact (QPC) in the presence 
of the Rashba spin-orbit interaction and a parallel mag¬ 
netic held. The quantum states near the center of the 
constriction are described by a one-dimensional massive 
Dirac equation, with the Rashba constant and magnetic 
held playing the roles of the speed of light and Dirac 
mass respectively. The emergent fermion and antifermion 
states behave differently in the channel, with the latter 
falling into the center of the constriction, forming reso¬ 
nant states. This process may be considered a manifesta¬ 
tion of Schwinger pair production, an effect which is well 
known in high-energy physics^^^— and has recently been 
investigated in the context of emergent two-dimensional 
Dirac system o^^i^^ . When the applied magnetic held is 
strong, the resonances become bound and generate a con¬ 
centration of spin current in the channel. The properties 
of these states are highly sensitive to the strength of the 
spin-orbit interaction and the applied magnetic held, im¬ 
plying that they can be controlled in situ by modulation 
of external helds. These results suggest an experiment 
in which repeated capture and release of particles inside 
the channel leads to a net spin polarization in the leads, 
transforming the channel into a spin pump and an in¬ 
jector of spin current which operates in the absence of a 
source-drain bias. While this effect exists in principle for 
both n and p type systems, consideration of experimental 
parameters for existing systems suggests that hole quan¬ 
tum point contacts provide an ideal candidate for the 
realization of this effect. 

A particle (either an electron or hole) moving ballisti- 
cally through a QPC in the presence of a Rashba spin- 
orbit interaction and parallel magnetic field is de¬ 
scribed by scattering states of the Hamiltonian 

= ^ - aPxOy - Pax + U(x) (1) 


where m is the band mass, a is the Rashba coefficient, 
V (cc) is a smooth scattering potential with maximum C/q 
at the center of the QPC [x = 0) and P = ^gpsBx- The 
spin-orbit interaction generates, in combination with the 
applied field, a dispersion consisting of spin-split bands 
given by 

^2 j ^2 _ 

e± = ^ ± + /32 . (2) 

For magnetic fields below a critical value P < Pc = ma^, 
the lower band has a “mexican hat” shape, with a local 
maximum at fc = 0; the dispersion in this case is shown 
in Fig. la. Note that the shape of the upper and lower 
branches near fc = 0 is governed by an anticrossing cre¬ 
ated by the combination of the Rashba interaction and 
a Zeeman interaction 0 < P < Pc- Inside the constric¬ 
tion, the kinetic energy is suppressed due to the barrier 
potential U{x). For states with certain incident energies, 
the scattering structure of the wavefunction is dominated 
by dynamics near the anticrossing in a region near the 
top of the barrier. In this region the linear-in-momentum 
term in o becomes dominant; to leading order in px we 
obtain the Dirac equation 

{-apxay - Pax + U{x))ii{x) = Eip{x) . (3) 

We may identify the upper and lower branches of the dis¬ 
persion, e^, near the anticrossing with the positive en¬ 
ergy and negative energy Dirac states respectively. Note 
that the parallel magnetic field creates a tunable Dirac 
gap of energy 2/3. 

Transmission through the constriction differs qualita¬ 
tively in the situations P > Pc and P < Pc- In the former, 
the transmission probability in each band simply drops to 
zero as the Fermi energy successively crosses the bottom 
of each band, leading to the well-known double-step in 
the channel conductance which is the experimental signa¬ 
ture of Zeeman split subbands^i. The second situation, 
P < Pc is highly non-trivial due to the possibility of tun¬ 
nelling from a positive to a negative energy state when 
the barrier is sufficiently high, Uq > E -\- p. Tunnelling 
occurs at the sides of the potential barrier (indicated by 
dashed lines in Fig. lb) in the regions where the kinetic 
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FIG. 1. (Color online) Confinement of a Dirac fermion in a 
QPC, which occurs when the parallel magnetic field is below 
the critical value /3 = < Pc- a) The ID disper¬ 

sion near the anticrossing, with the bold lines showing the 
region in which Dirac fermion behaviour is expected to ap¬ 
pear. The two spin branches e'p are identified as the pos¬ 
itive and negative energy Dirac states respectively. Arrows 
indicate the direction of spin polarization, b) The potential 
barrier U(x) (lower curve), and the situation in momentum 
space (three upper figures). There are three regions which 
form the scattering wavefunction. In the asymptotic region 
I a: I > b, the kinetic energy lies in the upper (positive energy) 
branch. Inside the barrier, |a;| < a, the kinetic energy lies in 
the lower (negative energy) branch. Quantization of motion 
in the region bounded by vertical lines gives rise to a subbar¬ 
rier resonance. The energy is indicated by the horizontal line. 
The tunnelling regions are shown with dashed lines. 


energy lies within the Dirac gap. The result is the ap¬ 
pearance of a particle in a negative energy state at the 

top of the barrier, which is characterized by a negative 

,2 

effective mass, < 0 (or equivalently, positive effec¬ 
tive mass but opposite electric charge) and therefore sees 
the repulsive potential as a quantum well. This demon¬ 
strates that there exist two types of charge carriers in the 
channel, which exhibit qualitatively different behaviour 
in the electrostatic potential U{x): the positive energy 
states in the branch fall away from the barrier, while 
the states in the antifermion branch ef fall towards the 
barrier. This process is illustrated in Fig. lb. The lower 
curve shows the potential U{x), with the tunnelling re¬ 
gion shown with dashed lines. The scattering wavefunc¬ 
tion is formed by two pairs of turning points x = ±6 
(assuming a symmetric barrier) and x = ±a. The situa¬ 
tion in momentum space is shown in the three hgures in 
the upper part of the same panel. In the regions |a;| > b, 
the kinetic energy (shown with a horizontal line) lies in 
the upper branch, and the particle accelerates away from 
the barrier. In the region \x\ < a, the kinetic energy lies 
in the lower branch. In this region the particle behaves 
as an antifermion; it falls into the barrier becomes con¬ 
fined. Quantization of motion between the turning points 
X = ±a gives rise to a resonant state; the energy of the 
resonant state is shown by the horizontal line intersecting 
the lower curve in Fig. lb. Note that the kinetic energy 


of the “antifermion” must remain above the bottom of 
the lower branch Emin in order for free motion to per¬ 
sist in the confinement region; the particle is expected 
to behave as a Dirac fermion when E — U{x) > Emin- 
The range of energies in which resonances are expected 
to exist is given by 

E + \Emin\ > Uq > E + P , Ennin =-^- 2ma^ ' 

""“( 4 ) 


The existence of two species of carriers which acceler¬ 
ate in opposite directions in the same electric field may be 
regarded as a manifestation of Schwinger’s mechanism for 
the spontaneous generation of fermion-antifermion pairs 
in strong electric fields^. While ubiquitous in quantum 
electrodynamics22i2i, direct observation of pair produc¬ 
tion requires electric fields eE > > 10^®V cm“^ 

which to date have only been accessible in high-energy 
collisions'^. In the condensed matter context, this ef¬ 
fect has been studied theoretically in graphen o^^i^^ and 
is closely related to the phenomenon of atomic col¬ 
lapse recently observed in the material^. In our case, 
the Schwinger mechanism distinguishes the localized 
“Dirac” state from quasibound states arising from dou¬ 
ble barriers, edge or multichannel effects as previously 
studie d^^d^~ — , and as we shall see later, gives them prop¬ 
erties which will enable them to facilitate pumping of spin 
across the channel. 

The transmission probability for a state incident in the 
upper branch may be derived by consideration of the 
explicit structure of the scattering wavefunction. The 
wavefunction consists of an asymptotic wave with pos¬ 
itive energy which undergoes reflection at the sides of 
the barrier, x = ±6, and a negative energy part cor¬ 
responding to free motion inside the barrier, |a;| < a. 
The positive energy and negative energy components of 
the wavefunction are coupled by tunnelling in the regions 
a < I a; I < 6, leading to a finite lifetime which is propor¬ 
tional to the inverse tunnelling rate. Approximating the 
barrier by a linear function, U'{x) = A, the resonant 
width is given by oc where the exponent 7 is 

given by Schwinger’s formulapi^ upon identification of /3 
with the rest energy and a with the effective speed of 
light: 


2aXh 


( 5 ) 


Accounting for periodic motion in the confinement re¬ 
gion |a;| < a and tunnelling, we obtain explicitly for the 
transmission probability 


T = 


g 27 r 7 _j_ (g 27 r 7 _ f kdx—iSip ' 


( 6 ) 


where j> kdx is the semiclassical phase acquired by the 
“antifermion” over one period of semiclassical motion. 
In the tunnelling regions, the solution of the Dirac equa¬ 
tion (|3]) may be expressed analytically in terms of the 













N = 0 


parabolic cylinder functions, which leads to a phase shift 
Sip = 2 Argr(i 7 ) — dy [ 1112^7 — l] + This phase shift 
appears in the WKB quantization condition for the ex¬ 
istence of a standing wave, j> kdx = 27r(n + 5 + Sip) 
at energies corresponding to the location of the Breit- 
Wigner resonances. We may regard m as the general 
expression for the transmission probability for a one¬ 
dimensional Dirac fermion in a smooth potential bar¬ 
rier. In the limit of strong magnetic fields /3 ^ V ah\, 
the resonant states become bound. In the opposite limit 
7 —>■ 0 , resonant tunnelling occurs over a broad range of 
energies, reflecting the Klein paradox for ultrarelativis- 
tic Dirac fermions^li^. Note that we have considered a 
state with asymptotic energy in the upper branch, ; in 
the semiclassical approximation a particle incident from 
the lower branch does not undergo reflection at the 
Dirac gap and the corresponding transmission probabil¬ 
ity is non-resonant. 

Approximating the QPC potential as a parabolic bar¬ 
rier, U{x) PS Uq— "^‘^2 ^ ; the resonant spectrum is equiv¬ 
alent (upon reversing the sign of the energy) to that of 
a harmonic oscillator with mass \m*\ and oscillator fre¬ 
quency U* = 

En = Uo- (3 - Uj*{n + ]^) . (7) 


The n = 0 mode possesses the highest energy, with higher 
modes forming an inverted tower of oscillator states ex¬ 
tending downward in energy. In deriving the spectrum 
© we assume that the lower band may be approximated 
by a quadratic dispersion e^, Ri —j3 -I- so that in 
the semiclassical picture the resonant states correspond 
to simple harmonic motion confined between two turning 
points X = ±a. (In this limit we also have Sip —>■ 0.) The 
spectrum terminates at finite n due to the condition (jd]) . 
This condition may be expressed in terms of parameters 

C = j 1 = as 

^ ma.‘^ ’ ' moc-^ 


C< 


I 


2,71 -|-1 


(1 — 77)27^2 


( 8 ) 


The number of quasibound states predicted by ([5]) in dif¬ 
ferent regions of the space of parameters (rj, () is shown 
in Fig. 2. In order to observe resonances we require 
the variation of the potential inside the channel to be 
smooth compared to the energy scale of the spin-orbit 
interaction, ^ <g; 1. The optimal regime, in which a large 
number of particles is trapped in the constriction, occurs 
when rj ~ Note that the number of quasibound states 
is much less sensitive to the size of the magnetic field 
than to the shape of the confining potential. 

It should be noted that number of states given by © 
may undergo modification in the presence of electron- 
electron interactions. In the Hartree approximation, the 
presence of a localized state at the center of the QPC 
both renormalizes the barrier height Uq and makes the 
barrier sharper near the top (increasing w), which would 
naively be expected to reduce the number of resonances. 
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FIG. 2. (Color online) Predicted number of quasibound states 
for parameters 77 , given by Eq. ®. 


At the same time, however, the spin-orbit coupling may 
be enhanced significantly via the exchange interaction in 
a similar manner to the exchange enhancement of the g- 
factor in quantum wire o^^’^° . Thus a detailed calculation 
beyond the Hartree approximation is required to account 
for the modification to © by interactions. 

Let us now address the possibility of observing reso¬ 
nant Dirac states in the typical experimental situation. 
Near pinchoff, the barrier height is tuned to the Fermi 
energy, Uo ~ Ep- We may parametrize the barrier in 
terms of the QPC length, Uq = ^muj^(^)^, to obtain 
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2 , fikF s2 
kpl ' ma 


(9) 


The value of I is limited by the ballistic mean free path, 
which is approximately one micron in 2D GaAs. Taking 
I = I/xm, at typical experimental density n = 10 ^^cm“^ 
in the 2D regions, we have kpl ~ 10^; in the strongly 
spin-orbit coupled materials InSb and InAs, the conduc¬ 
tion band effective mass (m ps 10~^me) is prohibitively 
small despite the strong Rashba effect, since ~ 10^ 
at typical electric fields a (x E^ = 10®Vm“^. While 
in principle the ID band-mass may be renormalized by 
inter-subband coupling, the initial value C ^ 10 ^ would 
require a renormalization by more than two orders of 
magnitude. 

Alternatively, we may consider p-type systems which 
possess a much larger value of m and a significantly 
stronger spin-orbit interaction^. In p-type point con¬ 
tacts, the size of the ID Rashba coefficient may be re¬ 
lated to the Rashba coefhcient in the 2D hole gas. The 
latter may be determined from the zero-held spin split¬ 
ting observed in magnetic oscillations; we obtain val¬ 
ues 0.6 < < 1 from such studies in GaAs-AlGaAs 

heterojunctions^i^. Taking typical values^ for the 
hole density n = 10 ^^cm“^ and m = 0.4me, we ob¬ 
tain 77 = -i, C = ^ at gBx = 3T for the lower value 
ma = 0.6hkp. In this regime © predicts N = 3. 














For comparison with the semiclassical solution de¬ 
scribed earlier, the scattering problem was also solved by 
explicit numerical integration of the Schrodinger equa¬ 
tion, with the corresponding QPC conductance shown in 
Fig. 3. The potential barrier was modelled as a Gaus- 
_ 

sian, U{x) = U^e ^ with w = 0.5/rm, and the con¬ 
ductance is shown for values of gB^ = IT, 2T, 3T (note 
that the typical 5 -factor in hole quantum point contacts 
is 5 ~ 0.5^). While the transmission problem was solved 
with the complete Hamiltonian o, the numerical solu¬ 
tion clearly exhibits the signatures of Dirac physics which 
were described earlier in semiclassical language based on 
the approximation ([3]). Results are shown in the case 
when the Fermi energy is high enough to intersect both 
branches of the dispersion in the asymptotic region; the 
conductance is given by the sum of transmission prob¬ 
abilities for particles injected in the upper and lower 
branches. 

For a particle injected in the lower branch ejT, the 
transmission probability smoothly decreases to zero when 
E — Uq < Emin, at which point there is reflection from 
the bottom of the lower band; the lower branch there¬ 
fore contributes a conductance G_ = ^ in the range (jH) . 

2 

The upper branch contributes a conductance G-\. = —T 
where the resonant transmission probability T is given 
by ([ 6 ]). Three resonances are shown as predicted by Eq. 
(| 8 |), with the roughly equal energy spacing reflecting the 
oscillator spectrum ©■ The lifetime of the rightmost 
resonance on the solid trace in Fig. 3 (indicated by the 
arrow) was calculated to be r = 7.4 x 10“^°s. The life¬ 
time, according to ([HI) and m, displays an exponential 
dependence on both a and /?, which may be used to pro¬ 
vide a highly sensitive measurement of the Rashba coef¬ 
ficient. Upon increasing the magnetic field from IT to 
3T, the resonances become significantly narrower (Fig. 
3). At 5T, the lifetime is ^ 10“®s, and at lOT, the life¬ 
time is ^ 10“^s. In this limit an additional broadening 
of the resonance is expected at finite temperature due to 
inelastic electron-electron collisions in the channel, nev¬ 
ertheless our analysis shows that at high magnetic fields, 
the resonance is effectively bound. Note that while vary¬ 
ing the magnetic field from IT to 3T significantly nar¬ 
rows the resonant peaks, it does not change the number 
of resonances. This distinguishing behaviour, which was 
derived earlier (see Eq. dS])) from the inverted harmonic 
oscillator spectrum 0 . permits a straightforward iden¬ 
tification of the effect in experiment. 



FIG. 3. (Color online) The QPC conductance in units of 

2 _ x2 

^ for a gaussian barrier U{x) = Uoe with parameters 
w = 50nm, Ef = O.BmeV, = 0.6, and values of magnetic 
field gBx = 3T (solid, red), 2T (dashed, blue), IT (dotted, 
green). The spin-resolved current density corresponding to 
the peak in the solid trace indicated by the arrow is shown in 
Fig. 4. 


presence of a Zeeman gap, only states of one chirality 
participate in the tunnelling process which leads to 
localization. In this sense the bound states considered 
here bear a strong similarity to chiral subgap (edge and 
surface) states in topological insulators^, although they 
are supported by a smoothly varying potential rather 
than an edge. The current densities of the spin-up and 
spin-down components of the wavefunction are shown in 
Fig. 4a for the wavefunction at the resonance indicated 
by the arrow in Fig. 3. 


It is straightforward to show that the localized spin 
current in our situation may be used to pump spin across 
the channel. Let us consider the spin states in the posi¬ 
tive and negative energy bands near the anticrossing. 
In the positive energy branch, states with positive mo¬ 
menta have spin tilted towards the —y axis, while states 
with negative momenta have spin tilted towards the +y 
axis. In the negative energy branch, the situation is op¬ 
posite (see Fig. la). Let us consider the component of 
spin current J^y which is given by 


T _ 1 rPx ", _ Px 01 


Let us finally consider the spin structure of the 
bound state. At low values of p (i.e. /3 <C ma^), 

the wavefunction inside the barrier is proportional to 
tp{x) oc where |±)y are spinors 

with polarization along the ±y axis. In such a state, 
both the current and magnetization are zero. Instead 
of a total spin, the state carries a total spin current 
Jxy.{x) = ^{vx,Sfj,} which is concentrated at the top 
of the barrier. The existence of a localized region of 
spin current is a consequence of the fact that, in the 


When a > 0 the first term is negative for a positive en¬ 
ergy state and positive for a negative energy state. Thus 
during the capture of a particle from the leads, the spin 
current of the particle is increased. In general, the bound 
state may be controlled by modulation of the parame¬ 
ters a, (3 which implies that the QPC can function as a 
spin transistor. We may perform, for example, a gedaken 
experiment in which the external magnetic field is first 
switched on slowly, resulting in the capture of a particle 
into the channel (assuming that the system is tuned so 
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FIG. 4. (Color online) Spin current and spin pumping via 
the bound state, a) The current densities (arb. units) in 
the spin-up and spin-down components of the wavefunction 
at resonance, Jx+ = 'ip+Jx'tp+ (red, positive curve) and Jx- = 
iplJxi’- (blue, negative curve), where ■0 = 'tp+\+)y+'>jj-\—)y. 
The barrier is shown by the dashed line. Since the state is lo¬ 
calized, the total current everywhere is zero, however the state 
possesses a spin current Jxy = , b) Upon adiabat- 

ically switching on the magnetic field, a particle is captured 
inside the barrier from the left lead, initially in a spin-down 
state, c) When the magnetic field is switched off instanta¬ 
neously, the resonance decays, either returning the particle to 
the left lead in a spin-up state, or transferring it to the right 
lead in a spin-down state. The mirror process occurs for a 
particle injected from the right lead. This process increases 
the number of spin-up carriers in the left lead and the num¬ 
ber of spin-down carriers in the right lead, i.e. it pumps spin 
across the channel in the absence of DC current. 


that the Fermi energy coincides with a resonance). Dur¬ 
ing this process a positive energy state in the reservoirs 
is converted into a negative energy state in the constric¬ 
tion, and the spin current of the system is increased. Let 
us now instantaneously switch off the external magnetic 


field. Since the external field is required for trapping, the 
bound state decays and the particle is transferred into ei¬ 
ther the left or right reservoir. However, in the absence 
of the external magnetic field, the operator of spin cur¬ 
rent (lion commutes with the Hamiltonian ([T|). Thus the 
initial spin current generated by adiabatically switching 
on the magnetic field persists even in the absence of a 
localized state. It is straightforward to show that this 
process either transfers a particle from one reservoir to 
the other without changing its spin, or returns a par¬ 
ticle to its original reservoir with a spin flip, as shown 
in Fig. 4b,c. Cycling the applied field in this manner 
pumps spin from one reservoir to the other, generating a 
non-equilibrium polarization in the leads. 

These results demonstrate that, aside from interac¬ 
tions and multichannel effects, the single particle, single 
mode transmission problem is in itself highly nontrivial 
in ID spin-orbit coupled systems. For the confinement 
mechanism described, only a strong spin orbit interac¬ 
tion, a smoothly varying repulsive potential, and an ex¬ 
ternal magnetic field are required. Furthermore, these re¬ 
sults illustrate a surprising manifestation of Dirac physics 
in ballistic ID channels, which is clearly indicated by 
the presence of Schwinger pair production, yielding an¬ 
tifermion states which undergo confinement and exhibit 
properties which are highly suitable for applications in 
spintronics. In the situation described, the parameters 
governing the Dirac equation are highly tunable, in con¬ 
trast to the situation in both high-energy physics and 
other emergent Dirac systems such as graphene and the 
topological insulators. In our case, the magnetic field 
plays a role analogous to the Dirac mass, the spin-orbit 
interaction plays the role of the speed of light, and the 
electric field is provided by the smooth ID potential in¬ 
side the constriction. The similarity and tunability of 
these energy scales allows for remarkable control over the 
properties of the Dirac states and enhances the versatil¬ 
ity of the system for applications to spintronics. The au¬ 
thor wishes to acknowledge O.P. Sushkov, J. Ingham, M. 
Veldhorst, I.S. Terekhov and A.R. Hamilton for valuable 
discussions. 
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